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Abstract We have simulated the decay of thermal counterflow quantum turbu- 
lence from a statistically steady state at r=1.9[K], with the assumption that the 
normal fluid is at rest during the decay. The results are consistent with the pre- 
dictions of the Vinen equation (in essence the vortex line density (VLD) decays 
as t~ l ). For the statistically steady state, we determine the parameter c% which 
connects the curvature of the vortex lines and the mean separation of vortices. A 
formula connecting the parameter %2 of the Vinen equation with C2 is shown to 
agree with the results of the simulations. Disagreement with experiment is dis- 
cussed. 

PACS numbers: 67.25.dk,67.25.dg 
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1 Introduction 

Superfluid 4 He at a finite temperature behaves as an an intimate mixture of two 
£C) | components: an inviscid superfluid component (density p s and velocity v s ); and 

7— I ■ a viscous normal component (density p n and velocity v n ). Quantum turbulence 1,2 

takes the form of a tangle of quantized vortices in the superfluid component; the 
normal fluid may or may not be turbulent. Thermal counterflow is an internal 
J"""-. 1 convection in which the superfluid component flows towards a source of heat, 
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while the normal fluid, carrying all the entropy, flows away from it. There is a 
relative velocity v ns = V n — V s between the two fluids, with no net mass flow. When 
the counterflow velocity exceeds a certain critical value, laminar flow proves to be 
unstable and the superfluid component becomes turbulent. This was the first type 
of quantum turbulence to be identified experimentally in the late 1950s 3 . 

The vortices in quantum turbulence interact with the normal fluid, giving rise 
a force of "mutual friction" between the two fluids. The study of the character- 
istics of this mutual friction provided the experimental evidence for counterflow 
turbulence. The turbulence in the superfluid component can be characterized by 
the vortex line density (VLD), L(t), and it was suggested that the time-evolution 
of this density in thermal counterflow might be described by the (Vinen) equation 

f-Wi" 1 -*^. (i) 

where a and %2 are temperature-dependent parameters and K is the quantum of 
circulation. The first term represents the generation of counterflow turbulence and 
the second term represents its decay. This equation implies that counterflow turbu- 
lence is homogeneous, which appears indeed to be the case as long as the spacing 
between vortices is much smaller than the size of the channel in which the heat 
current flows. 

Since the 1950s, many experimental, theoretical and numerical studies of coun- 
terflow turbulence have been reported. Schwarz showed by numerical simulations 
with the vortex filament model 4 that a self-sustained statistically steady-state form 
of turbulence can be maintained in the superfluid component 5 by the mutual fric- 
tion, and a more satisfactory simulation based on a full Biot-Savart simulation has 
recently been performed by Adachi et al. 6 . In both simulations flow of the normal 
fluid is assumed to be laminar. The steady-state turbulence is realized when the 
generation and decay terms in Eq. (QJ balance, so that L = L s t ea dy = a v ns> where 
7 is a temperature-dependent parameter. 

Counterflow turbulence in the steady state seems now to be well understood, 
but its decay, after the heat current is switched off, is less well understood. Ac- 
cording to Eq.([T]l, with v ns = 0, this decay should obey the equation 

I K I 

where Lq is the initial VLD. However, it has been known for a long time that 
in practice the decay proceeds in a more complicated (non-monotonic) manner 3 , 
although there is some evidence that Eq. (0 is obeyed at small times. It seems now 
to be well established 7 that for large times the decay proceeds as f -3 / 2 , implying 
that at this stage there is coupled motion of the two fluids dominated by eddies on 
a scale equal to the channel width 8 ; i.e. the normal fluid is also turbulent. 

In this paper we report the results of simulations of the decay of counterflow 
turbulence, starting from a statistically steady state, in order to check the predic- 
tions of Eq. (QJ. We assume that the normal fluid remains at rest during the decay. 
The validity of our simulations is limited to the extent that our VLDs are somewhat 
smaller than those encountered in practice, larger VLDs being hard to handle, but 
within this limitation our results prove to be consistent with the predictions of Eq. 
(QJ, with the expected value of X2- 
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2 The simulation 

We use the vortex filament model 4 ' 5 and represent a vortex as a string of discrete 
points s (%,t), where t is time and S, is arc length. At K the velocity s'o of the 
filament at the point s is given by 



(3) 



where the prime denotes differentiation with respect to arc length l + and l- 
are the lengths of the two adjacent line elements after discretization, separated by 
the point s, and is the cutoff corresponding to the core radius. The first term 
shows the localized induction field arising from a curved line element acting on 
itself. The second term represents the non-local field obtained by carrying out the 
integral of the Biot-Savart law along the rest of the filament. The third term v s is 
an applied field. 

At finite temperatures, it is necessary to take into account of the mutual friction 
between the vortex core and the normal flow v n , so that the velocity s becomes 

s = s + as' x (v n - s ) - a's' x [s' x (v n - i )] , (4) 

where a and a' are temperature-dependent friction coefficients 4 , and So is calcu- 
lated from Eq. (|3). 

Some important quantities that are useful for characterizing the vortex tangle 
can introduced 5 . The VLD is 

L=^f d$, (5) 

where the integral is performed along all vortices in the sample volume Q.. In 
the experiments 3,7 , the VLD was obtained from the attenuation of second sound 
propagating perpendicular to the counterflow v ns . Because mutual friction acts 
only at right angles to a vortex line element, the directly measured quantity is not 
the real VLD given by Eq. ©, but rather a "measured VLD", given by 

Lp = hL sin2ed ^ (6) 

where 6 is the angle between a vortex line segment and and the direction r of 
propagation of the second sound. If the quantum turbulence is isotropic, L p = 
(2/3 )L. Anisotropy of the vortex tangle is also interesting and is described by the 
dimensionless parameters 



1 

" ~ i2L J / 
1 



\-(s'-f 



d$, (7) 
d%> (8) 



where hi and represent unit vectors parallel and perpendicular to the original 
heat flow (in our case h^ = z). Since r± is the same as the direction of propaga- 
tion of the second sound, then L p /L = I±. Symmetry generally yields the relation 
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/ii/2+ij_ = 1. If the vortex tangle is isotropic, the averages of these parameters 
are 7|| =7^ =2/3. At the other extreme, if the tangle consists entirely of curves 
lying in planes normal to v ns , I\\ = 1 and 7± = 1 /2. 

In this study, calculations are performed under the following conditions. The 
numerical space resolution is A^ = 8.0 x 10 _4 [cm], the time resolution is At = 
1.0 x 10~ 4 [s], temperature is 1.9[K], and the computing box size is 0.1 x 0.1 x 
0.1 [cm 3 ]. We use periodic conditions in all directions. 



3 Results 

3.1 Turbulence in the statistically steady state 

Our decays start from a statistically steady state, and it is instructive first to de- 
scribe the character of this state and how it is generated 6 . 

First, we place in the box six vortex rings moving towards its centre, and we 
impose the relative velocity, v ns = 0.52[cm/s], along the z-axis at T=1.9[K]. The 
vortices make many reconnections and the VLD increases until it becomes statis- 
tically steady. The VLD of this steady state turbulence is Lq ~ 6.04 x 10 3 [1 /cm 2 ] 
and it satisfies the equation L = 'f'V^. 

We investigated the relationship 5 between the curvature of the vortex lines and 
the mean distance between vortices, 1 = 1 /L l l 2 , 

where the left hand side is the mean square curvature of the vortex lines and C2 
is a temperature-dependent parameter. The parameter C2 tends to decrease with 
increasing temperature, because the high velocities associated with small radii of 
curvature are more strongly damped by mutual friction at high temperature. 
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Fig. 1 (Color online) (a) Dependence of c-i on time in the statistically steady state. The parame- 
ter C2 fluctuates around 1.99. (b) Distribution of c%. The peak value is cj = 1 .96 and the standard 
deviation is a Cj = 0.38. 
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Figure Q] (a) shows that the parameter C2 is statistically constant and the time 
average of C2 is 1.99. Apparently, the fluctuation of C2 comes from the fluctuations 
in the statistically steady state of the VLD and the mean square curvature. We plot 
the probability density distribution of C2 in Fig.[T](b) obtaining C2 = 1.99±0.38, 
where the error is a standard deviation. The dependence of C2 on temperature will 
be reported in another paper. 



3.2 Decay of the counterflow turbulence 

After we obtained the steady state, we turned off the relative velocity, so that 
the turbulence can decay. It is necessary to study the statistics of the decaying 
turbulence, because the relatively small initial VLDs lead to fluctuations of the 
decaying VLD that are not negligible in one decay. Therefore we simulated 41 
decays from slightly different initial conditions, and then averaged the results. 
Figure [3] shows these averaged results, in the form of plots of L and l/L against 
time. Red circle points are the real VLDs and green triangle points show how the 
"measured VLDs" are expected to behave. 

Figure [2] shows the configuration of vortices in the decaying turbulence at 
various times. As we demonstrate below, the configuration is initially anisotropic 
(a), but after about 0.5 [s] it becomes isotropic (b), until at times greater than about 
3[s] the line density becomes so small (line spacing comparable with the size of 
the box) that the decay is strongly influenced by the walls of the box. 




Fig. 2 (Color online) Configuration of the vortex tangle at (a) J=0.0[s], (b) t ~0.5[s], and (c) 
/ >3.0[s]. 



We see from Fig.|3]that an equation of the form of Eq. (0 is obeyed as long as 
the VLD is not so small that the decay is influenced by the walls of the box. This 
result applies to both the real VLD and the measured VLD. 

It was argued in reference 9 from a consideration of the total rate of energy 
dissipation due to mutual friction that %2 an d C2 can be expected to be related by 
the equation 

X2=°f^, (10) 
2 go 

where Ro = l/(Lo c 2) an< ^ ls an initial VLD. In this way we can predict that 
at 1.9[K] %2 — 5.48 ±0.88. We see from Fig.[3]that our simulations support this 
prediction. 
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Fig. 3 (Color online) The decay of VLD (red circle points) and the measured VLD (green 
triangle points). These figures are averaged data. Figure (a) is a log-log plot of L against t, and 
Fig. (b) is a normal plot of 1 /L against t. The blue dotted line is the solution of Vinen equation 
(Xi = 5.48). The black dot-dashed line shows t power law. 



The way in which the degree of anisotropy changes during the decay is shown 
in Fig. [4] in which we have plotted the ratio L p /L and the parameter /i| against 
time. For isotropic turbulence L p /L = Ih = 2/3. 
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Fig. 4 (Color online) The ratio of L p to L (red circle points) and the anisotropy parameter, /y 
(blue triangle points), plotted against time. The black dotted lines shows 2/3. 



As we have ready noted, the decay at times greater than ~3[s] is strongly 
influenced by the walls of the box, and this can be seen in the behavior of L p /L 
in Fig. [4] (a). Figure [4] (b) shows that the decaying turbulence remains strongly 
anisotropic until f ^0.5 [s]. In the interval between /=0.5[s] and ?=2.5[s], the av- 
erages of Lp/L and /|| are 0.65 ± 0.04 and 0.70 ± 0.08, both values being close to 
2/3, indicating that the turbulence has become approximately isotopic. 
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4 Conclusions 

We have carried out computer simulations of the decay of quantum turbulence, 
using the vortex filament model 4 '^, with the assumption that the normal fluid is 
stationary during the decay. The simulations start from a statistically steady state 
of thermal counterflow turbulence 6 (Fig.[2](a)), generated by a flow in which v s — 
v n = 0.52[cm/s] at a temperature of 1.9[K]. The results prove to be consistent with 
the predictions of the Vinen equation (essentially the vortex line density decays 
with time as f _1 in the limit of large times, as long as the line spacing is small 
compared with the box size). A parameter C2, which describes the relationship 
between the curvature of the vortex lines and the mean distance between vortices, 
was evaluated for the steady state, and it is shown that the value of a parameter 
%2 that appears in the Vinen equation is correctly predicted in terms of C2 by a 
formula derived by Vinen and Niemela 9 . The turbulence is markedly anisotropic 
in the steady state counterflow, but it becomes isotropic as the decay proceeds. 
The validity of these results at temperatures other than 1.9[K] will be investigated 
in future work. 

The type of decay that emerges from our simulations does not agree with ex- 
periment, as has been suspected for many years. The origin of the discrepancy has 
been discussed by Schwarz and Rozen 10 and by Barenghi et al.—. Our results do 
not add support to the idea of Barenghi et al. that the discrepancy is associated 
in part with the transition during decay from an anisotropic to an isotropic vortex 
configuration. We suspect that the discrepancy is due primarily to the existence 
of coupled motion of the two fluids during the decay (perhaps even in the steady 
state), as suggested by Schwarz and Rozen. 
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